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imization, which takes fractional-orders and parameters as its particular in- 
dependent variables. And the objective is to find optimal combinations of 
fractional-orders and systematic parameters by DE in the predefined inter- 
vals for fractional order chaotic systems such that the objective function is 
minimized. Simulations are done to estimate a series of non-commensurate 
and hyper fractional chaotic systems. The experiments' results show that 
the proposed inversion mechanism for fractional-order chaotic systems is a 
successful methods with the advantages of high precision and robustness. 
Keywords: Non-commensurate and hyper fractional-order chaotic system, 
general mathematical model, inversion mechanism, unknown fractional 
orders and systematic parameters, non-negative special functions' 
minimization 
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1. Introduction 

Recently, the applications of fractional differential equations, whose non- 
linear dynamics are described by a powerful tool with the concept of frac- 
tional calculus [1-9], began to appeal to related scientists [10-29] in following 
areas, bifurcation, hyperchaos, proper and improper fractional-order chaos 
systems and chaos synchronization [8, 12-14, 16-18, 24-32]. 

For the fractional-order chaos systems, most of the control and synchro- 
nization methods are generalized conclusive but much more complicated ap- 
proaches from the methods for normal and hyper chaos system. However, 
when considering the controlling and synchronization for the fractional-order 
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chaos systems, there are some systematic parameters and fractional orders 
are unknown. It is difficult to identify the parameters in the fractional-order 
chaotic systems with unknown parameters [16, 17, 33, 34]. 

The process to get the exact values of uncertain orders and parameters 
for the fractional order chaotic systems is called system inversion mecha- 
nism. Hitherto, there have been some approaches in system inversion for 
fractional-order chaos systems. For instances, synchronization for fractional- 
order chaos[33] and fractional order complex networks[34]. However, the 
design of controller and the updating law of parameter identification is still 
a hard task with techniques and sensitivities depending on the considered 
systems. And the non-classical way via artificial intelligence methods, for 
examples, differential evolution [16] and particle swarm optimization [17]. In 
which only the commensurate fractional order chaos systems and simplest 
case with one unknown fractional order for fractional-order chaos systems 
are discussed. And there exist basic hypotheses in traditional non-Lyapunov 
estimation methods[10, 16, 17]. That is, the parameters and fractional or- 
ders are partially known or the known data series coincide with definite 
fractional equationas forms / = (/i, f2, f n ) of fractional chaotic differ- 
ential equations except some uncertain parameters and fractional orders 
6 = # 2 , 8 n , qi, (72, q n ). In which, the non- commensurate cases with 
different fractional orders {qi, % = l,...,n} are not included in the above 
non-Lyapunov ideas or not fully discussed either. 

However, to the best of our knowledge, few work in non-classical way has 
been done to the parameters and fractional orders' inversion estimation of 
non-commensurate and hyper fractional-order chaos systems [35] so far. 
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The objective of this work is to present a novel simple but effective in- 
version mechanism of uncertain fractional-orders and parameters based on 
differential evolution algorithms (DE) to estimate the non-commensurate and 
hyper fractional order chaotic systems in a non-Lyapunov way. And the illus- 
trative inversion simulations in different chaos systems system are discussed 
respectively. 

The rest is organized as follows. Section 2 give a simple review on non- 
Lyapunove parameters inversion mechanism for fractional-order chaos sys- 
tems. In Section 3, a general mathematical model for fractional chaos pa- 
rameters inversion mechanism in non-Lyapunov way are analyzed. In section 
4, firstly the DE are introduced briefly. Then a novel methods with proposed 
united model based on DE is proposed to estimate the fractional chaos sys- 
tems. And simulations are done to a series of different non-commensurate 
and hyper fractional order chaotic systems. Conclusions are summarized 
briefly in Section 5. 

2. Non-Lyapunove parameters' inversion mechanism for fractional 
order and normal chaos systems 

We consider the following fractional-order chaos system. 

' a D q tyi (t) = h(yi(t),-,yn(t),Y (t),8); 

a D q t y 2 (t) = f 2 ( yi (t),...,y n (t),Y (t),e); 

< (1) 
a D q t y n (t) = f n ( yi (t),...,y n (t), Y (t), 6) . 

L = (2/1,2/2, --,2/n) 

Let L{t) = (yi(t),y 2 (t), ...,y n (t)) T G 5i n denotes the state vector. 9 = 

4 



9 2 , 9 n ) T denotes the original parameters, q = (q 1} q 2 , q n ), (0 < % < 
l,i = 1, 2, n) is the fractional derivative orders. 

Normally the function / = (fx, f 2 , /„) is known. And the 9 = (9\, 9 2 , 9 n ) 
q = (q±, q 2 , q n ) will be the parameters to be estimated. Then a correspon- 
dent system are constructed as following. 



a DWi (t) = fi (yi (t), y n {t), Y (t),9) ; 

a D q t y 2 (t) = h (Vi(t), Ut), Y (t),6) ; 

a D q t y n (t) = f n y n {t), Y (t), 9) 
L = (yi,y 2 , ...,y n ) 



(2) 



where L(t), 9, q are the correspondent variables to those in equation (1), and 
function / are the same. The two systems (1) (2) have the same initial 
condition Y Q (t). 

Then the objective is obtained as following, 



q)* = arg min F = arg min 2_. \\L (t) — L (t) 

(&,q) (0,q) t 11 



(3) 



When some the fractional chaotic differential equations / = (fi, f 2 , f n ) 
are unknown, the objective will be, 



(fi,f2,--,fn)* = arg min F 

(/l,/2,...,/n) 



(4) 



Equation (4) is fractional-order chaos' inversion mechanism called recon- 
struction. In Reference [36] a novel method was proposed to reconstruct the 
unknown equations (/i, f 2 , f n ) based on an united mathematical model. 
However, for the united mathematical model [36], to be identified is only 
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{hi f2, fn) instead of q. That is, the left part a DlY{t) of the equation (2) 
are not included. 

Therefore, to estimate the q of the equation (2) with unknown system- 
atic parameters 9 is still a question to be solved for parameters and orders 
estimation of non-commensurate and hyper fractional-order chaos systems. 

Now we take the non-commensurate fractional order Lii system (5) [22, 37] 
for instance. 

' D?x(t) = a(y(t)-x(t)); 

oDfy (t) = -x (t) z (t) + cy (t) ; (5) 
^ vDfz(t)=x(t)y(t)-bz(t). 
when (a, b, c) = (36, 3, 20), (q u q 2 , q 3 ) = (0.985, 0.99, 0.98), initial point (0.2, 0.5, 
Lii system (5) is chaotic. When (a, b, c), (q±, q%, q^) is unknown, we have to 
estimate them. 

Secondly the objective function is chosen as: 

T-h 

p = F(a,b,c, qi ,q 2 ,q 3 )= ^ ||y (t) — Y (t)||^ (6) 

t=0-h 

The objective function for system (5) is shown as Figure 1. 

Then the problems of inversion for chaotic system are transformed into 
that of nonlinear function optimization (6). And the smaller p is, the better 
combinations of parameter (a, b, c, q±, q 2 , (fe) is. The independent variables of 
these functions are 9 = (a, b, c, qi, q 2 , (fe). 

For the normal chaos system, parameters inversion of chaotic system is 
the chief task to be resolved and of vital significance, especially for some 
new hyperchaotic systems recently proposed [38-46]. And these systems can 
be thought as special case of the above fractional order chaos systems. And 
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Figure 1: The objective function of fractional order Lii system 

series of estimation methods are proposed based on genetic algorithms (G A), 
chaotic ant swarm algorithms, particle swarm optimizations(PSO), quantum 
PSO, artificial bee colony algorithms etc[47-61]. 

3. General mathematical model for fractional chaos inversion in 
non-Lyapunov way 

In this section, a general mathematical model for fractional chaos param- 
eters identification in non-Lyapunov way is proposed. 

The general forms of fractional order chaos systems is (1). We consider 
that (/i, / 2 , f n ) is to be reconstructed and (qi, q2, q n ) is to be identified 
in (2). 

To have simple forms, we take a = 0. The continuous integro-differential 
operator[37, 62] is defined. And for the continuous function /(£), the nu- 
merical solution method are used here [37, 62], which was obtained by the 
relationship derived from the G-L definition to resolve system. 
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Then the novel objective function (fitness) equation (7) in this paper come 
into being as below. 

T-h 

F =H\\ l - L \\ 2 ( 7 ) 

t=0h 

It should be noticed here that the independent variables in function (7) 
are not the parameters as in (6) but the special variables, for instance, as 
functions (/i, f 2 , f n ) and fractional orders (qi, q 2 , ■■, q n ). That is 

((/i,/2, -,/n), (?i,?2, -,?n)) = _ _ arg min F (8) 

((/l>/s,-./n)i(?l,«a,".,Si)J 

Equation (8) is the crucial turning point that changing from the parame- 
ters inversion into fractional equations and orders inversion, in other words, 
both fractional order estimation and fractional chaos systems' reconstruction. 

Now, with the definitions of equations (7) and (8), a novel general mathe- 
matical model for fractional chaos parameters identification in non-Lyapunov 
way is coming into being. 

It can be concluded that the parameters' estimation of fractional order 
chaos system[16, 17] is a special case of fractional order chaos reconstruction 
here as (8). 

However, it should be emphasized here that, it is not easy to reconstruct 
the fractional order differential equations and identify the fractional orders 
together. And only the simplest case that with definite q = q x = q 2 = ... = q n 
are discussed [36]. 

And for a special cases of the general model (8) that (qi ^ q 2 ^ ... ^ 
q n ) and systematic parameters 9 2 , 9 n ) G f 2 , / n ) are unknown 
for non-commensurate and hyper fractional order chaos system, to our best 
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of knowledge, no such inversion mechanisms have been done. Thus, it is 
necessary to resolve the following (9) in non-Lyapunov way. 

(q,9)* = argminF (9) 

(1,8) 

4. A novel inversion mechanism based on differential evolution al- 
gorithm 

The task of this section is to find a simple but effective inversion ap- 
proach for unknown q and systematic parameters in equation (9) of for non- 
commensurate and hyper fractional-order chaos based on Differential Evolu- 
tion (DE) algorithm in non-Lyapunov way. 

4-1. The Main Concept of Differential Evolution Algorithm 

Differential Evolution (DE) algorithm grew out of Price's attempts to 
solve the Chebychev Polynomial fitting Problem that had been posed to 
him by Storn [63]. A breakthrough happened, when Ken came up with 
the idea of using vector differences for perturbing the vector population. 
Since this seminal idea, a lively discussion between Ken and Rainer and 
endless ruminations and computer simulations on both parts yielded many 
substantial improvements which make DE the versatile and robust tool it is 
today [63-66]. 

DE utilizes M n-dimensional vectors, Xj = (xn, • • • , Xj n ) G S, i = 1, • • • , M, 
as a population for each iteration, called a generation, of the algorithm. For 
each vector x\ G ^ = (X^, X^\ ■ ■ ■ , xf^) : i = 1, 2, • ■ ■ , M, there are three 
main genetic operator acting. 
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For each individual, to apply the mutation operator, firstly random choose 
four mutually different individual in the current population ,Xr 2 ,Xr 3 {r\ ^ 
r 2 7^ ^3 7^ i) -Then combines it with the current best individual X^ est to get 
a perturbed vector V = (V\, V 2 , ■ ■ ■ , V n ) [63, 67] as below: 



V 



Xlf ] + 0.5(CF + 1) • (x r (G) + X r (G) - 2Xi G) ) , if rand{0, 1) < 0.5 

X4 G) + CF ■ (x^ ] - X r ( 2 G) ) , otherwise 

(10) 

where CF > is a user-defined real parameter, called mutation constant, 
which controls the amplification of the difference between two individuals to 
avoid search stagnation. 

(G) 

Following the crossover phase, the crossover operator is applied on X- . 
Then a trial vector U = (Ui, U2, ■ ■ ■ , U n ) is generated by: 

V m , if (rand(0, 1) < CR) or (m = k), 

Um — \ (fA (11) 

X\ G 2 , if {rand{0, 1) > CR) and {m^k). 
in the current population [63], where m — 1, 2, • • ■ , n, the index k e {1, 2, • ■ ■ , n} 
is randomly chosen, CR is a user-defined crossover constant [63, 67] in the 
range [0, 1]. In other words, the trial vector consists of some of the compo- 
nents of the mutant vector, and at least one of the components of a randomly 
selected individual of the population. 

Then it comes to the replacement phase. To maintain the population 

(G) 

size, we have to compare the fitness of U and X\ , then choose the better: 

X (c + i )=< ! U,if F(U) < F(XP), (i2) 
X„- G ^ , otherwise. 
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4-2. A novel unknown parameters and fractional orders inversion mechanism 

Now we can propose a novel approach for hyper, proper and improper 
fractional chaos systems. The pseudo-code of the proposed reconstruction is 
given below. 

Algorithm 1 A novel inversion mechanism based on differential evolution 

algorithms 
1: Basic parameters' setting for DE. 

2: Initialize Generate the initial population. 

3: while Termination condition is not satisfied do 

4: Evaluation Evaluate the fitness with Eq. (9) and remain the best 

individual. 
5: Mutation As in equation (10). 
6: Crossover As in equation (11). 
7: Replacement As in equation (12). 

8: Boundary constraints For each x^ G X t , k = 1,2, ...,D, if Xn is 
beyond the boundary, it is replaced by a random number in the boundary. 
9: end while 

10: Output Global optimum xsest 



4-3. Non- commensurate and hyper fractional order chaos systems 

To test the Algorithm 1, some different well known and widely used non- 
commensurate and hyper fractional order chaos systemsare choose as follow- 
ing. 

Example. 1. Fractional Lorenz system[13, 28, 68] for instance, which is 
generalized from the first canonical chaotic attractor found in 1963, Lorenz 
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system[69]. 



(13) 



a Dfx = a ■ (y - x); 
a D q t 2 y = 7 • x - x • z-y; 
ojyf'z = x ■ y — b ■ z. 
L = (x,y,z) 

where qi, q 2 , q% are the fractional orders. When q 2 , q^) = (0.985, 0.99, 0.99), 
a = 10, 7 = 28, b = 8/3, a = 0, intimal point (0.1,0.1,0.1), it is the non- 
commensurate chaotic system[68]. The objective function for system (13) is 
shown as following. 




(o)(b,Sj, Ji) = (0.985,0.99,0.99) 




ffl(»,7,&) = (10,28,8/3) 



Figure 2: The objective function of fractional order Lorenz system 

Example. 2. Fractional order Arneodo's System (14) [37, 70]. 
' D?x(t)=y(t); 

< ,DTy{t) = z{t)- (14) 
^ Df*z (t) = -ftz (t) - /3 2 y (t) - f3 3 z (t) + fax 3 (t) . 

when (ft, ft, ft) = (-5.5,3.5,0.8,-1.0), (gi,g 2 , ?3 ) = (0.97,0.97,0.96), 
initial point (—0.2,0.5,0.2), Arneodo's System (14) is chaotic. 
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Example. 3. Fractional order Duffmg's system (15) [37]. 
D?x(t) = y(t); 

(15) 

Q Dfy (t) =x(t)- x 3 (t) - ay (t) + S cos (cot) . 

when (a,b,c) = (0.15,0.3,1), {q 1 ,q 2 ) = (0.9,1), initial point (0.21,0.31), 
Duffings system (15) is chaotic. 

Example. 4. Fractional order Genesio-Tesi's System (16) [37, 71]. 
' D?x(t)=y(t); 

< ,Dfy{t) = z{t)- (16) 
^ Q D?z (t) = -f3 lX (t) - (3 2 y (t) - (3 3 z (t) + fax 2 (t) . 

when (Pi, p 2 , P 3 , = (1.1,1.1,0.45,1.0), (gi,g 2 ,g 3 ) = (1,1,0.95), initial 
point (—0.1,0.5,0.2), Genesio-Tesi's System (16) is chaotic. 
Example. 5. Fractional order financial System (17) [37, 72] with the exact 
form of the differential equation oDf 3 z = f 3 (x, y, z) are unknown. 

D?x (t) = z(t)+x [t) (y (t) - a) ; 
< D?y(t) = l-by(t)-T t {t); (17) 
^ Q Dfz(t) = -x{t)-cz(t). 

when (a, b, c) = (1, 0.1, 1), (qi, q 2 , q 3 ) = (1, 0.95, 0.99), initial point (2, —1,1), 

financial System (17) is chaotic. 

Example. 6. Fractional order Lii system (5). 

Example. 7. Improper fractional order Chen system (18) [37, 73, 74]. 
' D?x (t) = a(y(t)- x (t)); 

< Dry(t) = (d)x(t)-x(t)z(t) + cy(t); (18) 
^ D?z(t)=x(t)y(t)-bz(t). 
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And when when (a, b, c, d) = (35,3,28,-7), (gi,g 2 ,<?3) = (1,1.24,1.24), ini- 
tial point (3.123,1.145,2.453), Chen system (18) is an improper chaotic 
system[18]. 

Example. 8. Fractional order Rossler System (19) [37, 75]. 



D?x(t) = -{y{t) + z(t)); 
D q t 2 y (t) =x{t) + ay (t) ; 
Dfz(t) = b + z(t) (x (t)-c). 



(19) 



when {a,b,c) = (0.5,0.2, 10), (gi,g 2 ,g 3 ) = (0.9, 0.85, 0.95), initial point (0.5,1.5,0.1), 
Rossler System (19) is chaotic. 

Example. 9. Fractional order Chuas oscillator (20) [76]. 



where 



oDfx (t) = a{y{t)-x (t) + (t) -W(w)x (t)) 
D?y(t)=x(t)-y(t) + z(t); 
Q DTz{t) = -Py{t)- 1 z{t)- 
D q t 4 w(t)=x(t); 



W(w) 



(20) 



a : 


\w\ 




\w\ 



when (a, f3, 7, (, a, b) = (10, 13, 0.1, 1.5, 0.3, 0.8), (q u q 2 , q 3 , q 4 ) = (0.97, 0.97, 0.97, 0.97), 
initial point (0.8,0.05,0.007,0.6), Chua's oscillator (20) is chaotic. 
Example. 10. Hyper fractional order Lorenz System (21) [77] . 



Dfx{t) = a{y{t)-x{t))+w{t)- 
Df 2 y (t) = cx (t) -x(t)z it) -y(t); 
D q t 3 z{t)=x(t)y(t)-bz(t); 
Dpw(t) = -y(t)z(t)+ 1 w (t); 



(21) 
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when (a,b,c,d) = (10,8/3,28,-1), (q 1} q 2 , q 3 , g 4 ) = (0.96,0.96,0.96,0.96), 
initial point (0.5,0.6,1,2), Hyper fractional order Lorenz System (21) is 
chaotic. 

Example. 11. Hyper fractional order Lii System (22) [78]. 

Dfx(t) = a(y it) -x it)) + w{t)- 
D?y(t) = -x(t)z(t) + cy(t); 
Dfz(t)=x(t)y(t)-bz(t) ] 
D q t 4 w (t) = x{t)z (t) + dw (t) ; 
when (a, b, c, d) = (36, 3, 20, 1.3), (q u q 2 , q 3 , q*) = (0.98, 0.980.98, 0.98), initial 
point (1, 1, 1, 1), Hyper fractional order Lii System (22) is chaotic. 
Example. 12. Hyper fractional order Liu System (23) [79]. 

Dfx (t) = -ax (t) + by (t) z{t) + z (t) ; 
Dr y (t) = 2.5y (t)-x(t)z(t); 
D q t 3 z (t) =x{t)y (t) - cz (t) - 2w (t) ; 
Q Dfw (t) = -d ■ x (t) . 
when (a,b,c,d) = (10,1,4,0.25), (q u q 2 , g 3 , g 4 ) = (0.9,0.9,0.9,0.9), initial 
point (2.4, 2.2, 0.8, 0), Hyper fractional order Liu System (23) is chaotic. The 
objective function for system (23) is shown as following. 
Example. 13. Hyper fractional order Chen System (24) [80]. 

D?x(t) = -a(y (t)-x(t))+w(t); 
D g t 2 y (t) = dx (t) -x(t)z (t) + cy (t) ; 
D?z(t)=x(t)y(t)-bz (i); 
^ D q t 4 w (t) =y(t)z (t) + rw (t) . 
when (a,b,c,d) = (35,3,12,7,0.5), (q lt q 2 , q 3 , g 4 ) = (0.96,0.96,0.96,0.96), 
initial point (0.5, 0.6, 1, 2), Hyper fractional order Chen System (24) is chaotic. 



(23) 



(24) 
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Figure 3: The objective function of Hyper fractional order Liu System system 



Example. 14. Hyper fractional order Rossler System (25) [75]. 



D q t 2 y (t) =x{t) + ay (t) + w (t) 
D?z(t) = x(t)z(t) + b; 
Q Dfw (t) = -cz (t) + dw (t) . 



(25) 



when (a,b,c,d) = (0.32,3,0.5,0.05), (91,92,93,94) = (0.95,0.950.95,0.95), 
initial point (—15.5, 9.3, —4, 18.6), Hyper fractional order Rossler System (25) 
is chaotic. 

Example. 15. A four- wing fractional order system[81, 82] both incommen- 
surate and hyper chaotic. 



Dfx x 



Dfx 2 
Dfx* 
D q t A Xi 



axi - x 2 x 3 + x 4 , 
-bx 2 + x x x 3 , 

XiX 2 — CX 3 + X1X4, 

-x 2 , 



(26) 
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when (a,b,c) = (8,40,49), (q u q 2 , g 3 , g 4 ) = (1,0.950.9,0.85), initial point 
(1, —2,3, 1) [81] , system (26) is chaotic. 

4-4- Simulations 

For systems to be identified, the parameters of the proposed method are 
set as following. The parameters of the simulations are fixed: the size of the 
population was set equal to M = 40, generation is set to 500, the default 
values CF = 1, CR = 0.85; The times of function evaluations are 20040. 
Table 1 give the detail setting for each system. 

Table 2 shows the simulation results of above fractional order chaotic 
systems. 

The following figures give a illustration how the self growing evolution 
process works by DE Algorithm 1. In which, Figures 4,5 ,6,7 8,9, 10,11 show 
the simulation evolution results of above fractional order chaotic systems with 
optimization process of objective function's evolution and the parameters and 
orders uncertain of above fractional order chaotic systems. 

From the simulations results of reconstruction above fractional order 
chaos system, it can be concluded that the proposed method is efficient. 
And from above figures, it can be concluded that the estimated systems are 
self growing under the genetic operations of the proposed methods. 

To test the performance of the proposed method Algorithms 1 , some 
more simulations are done to the four-wing incommensurate hyper fractional 
order chaotic system (26) in following cases A,B,C,D. In these cases, each 
with only one condition is changed according to the original setting for system 
(26). And the simulation results are listed in Table 3. 
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Table 1: Detail parameters stetting for different systems 



F-0 systems Unknown Lower boundary Upper boundary Step No. of samples 



Lorenz 


(c-,7,6, 91,92,93) 


5,20,0.1,0.1,0.1,0.1 


15,30,10,1,1,1 


0.01 


100 


Arneodo 




-6,2,0.1,-1.5,0.1,0.1,0.1 


-5,5,1,-0.5,1,1,1 


0.005 


200 


Duffing 


(a,b,c, 91,92) 


0.1,0.1,0.1,0.1,0.5 


1,1,2,1,1.5 


0.0005 


500 


Genesio-Tesi 


(0i,02,0a,04.,qi,q2,qs) 


1,1,0.1,0.1,0.5,0.5,0.1 


2,2,1,1.5,1.5,1.5,1 


0.005 


200 


Financial 


(a,b, c, 91,92,93) 


0.5,0.01,0.5,0.5,0.1,0.1 


1.5,1,1.5,1.5,1,1 


0.005 


200 


Lii 


(a, b, c, 91,92,93) 


30,0.1,15,0.1,0.1,0.1 


40,10,25,1,1,1 


0.01 


100 


Improper Chen 


(a,6,c, d, 91,92,93) 


30,0.1,20,-10,0.5,1,1 


40,10,30,-0.1,2,2,2 


0.01 


100 


Rossler 


(a, b, c, 91,92,93) 


0.1,0.1,5,0.1,0.1,0.1 


1,1,15,1,1,1 


0.01 


100 


ChuaM 


(a,/3,7,C,a,M) 


5,10,0.1,0.1,0.3,0.1,0.1 


10,20,1,2,0.3,1,1 


0.01 


100 


Hyper Lorenz 


(a, 6, c, d, 9) 


5,0.1,20,-2,0.1 


15,5,30,-0.1,1 


0.01 


100 


Hyper Lii 


(a, b, c, d, 9) 


30,0.1,15,0.1,0.1 


40,5,25,5,1 


0.005 


200 


Hyper Liu 


(a, b, c, d, 9) 


5,0.5,1,0.1,0.1 


15,1.5,10,1,1 


0.005 


100 


Hyper Chen 


(a, b, c, d, 9) 


30,0.1,10,0.1,0.1,0.1 


40,5,20,10,1,1 


0.005 


200 


Hyper Rossler 


(a, b, c, d, 9) 


0.1,0.1,0.1,0.01,0.1 


1,5,1,1,1 


0.005 


200 


System (26) 


(a, b, c, 91,92,93,94) 


5,38,45,0.5,0.5,0.5,0.5 


10,45,50,1,1,1,1 


0.005 


200 



a For the hyper fractional order chaos system, let q = q\ = q% = q% = q±. 



Table 2: Simulation results for different fractional order chaos systems 



F-0 system StD Mean Min Max Success rate a 



Lorenz 


1 


2271e-05 


2. 


3801e-05 


3. 


6823e- 


•06 


5. 


4275e- 


■05 


100% 


Arneodo 


6, 


,1944e-07 


1. 


4863e-06 


1. 


6606e- 


•07 


1. 


3293e- 


•06 


100% 


Duffing 


2. 


4378e-10 


2. 


,4788e-10 


2. 


,867e-ll 


1. 


7728e- 


•09 


100% 


Genesio-Tesi 


1. 


76e-06 


3. 


7268e-06 


6. 


6136e- 


•07 


1. 


0572e- 


•05 


100% 


Financial 


2. 


3322e-07 


4. 


5927e-07 


1. 


1207e- 


•07 


1. 


1711e- 


•06 


100% 


Lii 


1. 


8797e-05 


2. 


9371e-05 


5. 


2309e- 


•06 


1. 


0775e- 


•04 


99% 


Improper Chen 


2. 


5763e-03 


1. 


1358e-03 


6. 


6422e- 


•04 


1. 


8223e- 


■02 


97% b 


Rossler 


8. 


4114e-08 


1. 


5306e-07 


2. 


4063e- 


•08 


1. 


5141e- 


•07 


100% 


ChuaM 


1. 


6998e-06 


2. 


7689e-06 


1. 


2085e- 


•07 


7. 


9197e- 


•06 


100% 


Hyper Lorenz 


7. 


2263e-08 


9. 


5764e-08 


1. 


4287e- 


•08 


3. 


8946e- 


•07 


100% 


Hyper Lii 


5. 


9079e-08 


1. 


037e-07 


2. 


2241e- 


•08 


2. 


9391e- 


•07 


100% 


Hyper Liu 


1. 


0783e-09 


1. 


6299e-09 


1. 


3241e- 


•010 


6. 


2280e- 


•09 


100% 


Hyper Chen 


1. 


1798e-05 


2. 


0742e-05 


6. 


1987e- 


•06 


6. 


0546e- 


•05 


100% 


Hyper Rossler 


6. 


120le-09 


6. 


9673e-09 


1. 


7082e- 


■09 


4. 


9026e- 


•08 


100% 


System (26) 


2. 


9764e-03 


5. 


,7449e-03 


1. 


7212e- 


•03 


1. 


.6418-02 


90% 6 



a Success means the the solution is less than le — 4 in 100 independent simulations. 



b Success means the the solution is less than le — 2 in 100 independent simulations. 
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Generations Generations 



Figure 5: Evolution process for fractional order Arneodo system 
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Figure 6: Evolution process for fractional order Financial system 



100 200 300 400 500 

Generations 




200 300 
Generations 



Figure 7: Evolution process for fractional order Lii system 
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6, 




100 200 300 400 500 100 200 300 400 



Generations Generations 



Figure 8: Evolution process for fractional order improper Chen system 
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Figure 9: Evolution process for fractional order hyper Lii system 
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Figure 10: Evolution process for fractional order hyper Liu system 




Figure 11: Evolution process for fractional order hyper Rdsslcr system 



23 



• Case A. Enhancing the defined intervals of the unknown parameters 
and orders to [0, 10] x [30, 40] x [40, 50] x [0.1, 1] x [0.1, 1] x [0.1, 1] x [0.1, 1] 



• Case B. Minimizing the number of samples for computing system (26) 
from 200 to 100. 

• Case C. Changing the iteration numbers of Algorithms 1 from 500 to 
800. 

• Case D. Changing the population size of Algorithms 1 from 40 to 80. 



Table 3: Simulation results for system (26) 



system (26) 


StD 


Mean 


Min 


Max 


Success rate a 


NEOF 6 


Case A. 


3.6259e+02 


3.5359e+02 


4.3690e-03 


7.2477e+02 


20% 


20040 


Case B. 


5.2549e-05 


1.1430e-04 


3.8482e-05 


2.7105e-04 


100% 


20040 


Case C. 


5.8123e-06 


1.0498e-05 


1.7082e-06 


3.0149e-05 


100% 


32040 


Case D. 


3.3265e-03 


9.6975e-03 


2.9053e-03 


2.1081e-02 


64% 


40080 



a Success means the the solution is less than le — 2 in 100 independent simulations. 
b No. of evaluation for objective function 



Figure 12 show the coresspondent simulation results for system (26). 

From results of the Table 2, 3 and Figure 12, we can conclude that mini- 
mizing the number of samples for computing the system (26) as case B, en- 
hancing the iteration numbers as case C, the population size of Algorithms 
1 as case D, will make the Algorithms 1 much more efficient and achieve 
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(a) Process of one simulation (b) Results of F for all simulations 

Figure 12: Simulation results for system (26) 

a much more higher precision. However if the defined intervals of the un- 
known parameters of system (26) are enhanced, then the results will go to 
the opposite way. That is the success rate is from 90% to 20% as case A. 

And considering the No. of evaluation for objective function, it is that 
minimizing the number of samples for computing the system (26) as case B 
is the best way to achieve higher efficiency and precision. 

5. Conclusions 

The inversion mechanism put consists of numerical optimization problem 
with unknown fractional order differential equations to identify the chaotic 
systems. Simulation results demonstrate the effectiveness and efficiency of 
the proposed methods with the mathematical model in Section 3. This is 
a novel Non-Lyaponov way for fractional order chaos' unknown parameters 
and orders. 

The performance of the proposed method is sensitive to a series factors, 
such as the initial point for each fractional order chaos system, sample in- 
terval, number of points, and length of intervals for the unknown orders and 
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parameters. Actually, these also lead to the candidate system divergent. And 
they are not predefined randomly. A good combination of these is not easy 
to get. Some mathematical formula to get a good combination not by so 
many simulations will be introduced in the future studies. 

It should be noticed too many points for evaluating fractional order chaos 
system the individual represents are not worth. Because the most time con- 
sumption parts in the whole proposed method are to resolve the candidate 
systems. Some of these system are easy to solve. However when it comes 
with the some individuals with bad combinations of parameters and orders , 
the methods to resolve the fractional order chaos systems in Section 3 might 
not converge as shown in the simulations for system (26). Then the whole 
proposed method Algorithms 1 might get into endless loops. To avoid the 
endless loops, we introduce a forced strategy to assign all the NAN and in- 
finite numbers in the output as zero. Because the objective function (9) to 
be optimized is bigger than 0, so this forced strategy for assignment is rea- 
sonable. To achieve a fine balance between the performance of the proposed 
methods and having enough sample data for credibility, we take the number 
of the points as 100 — 200, according the existing simulations [47-61]. And 
the simulations in section 4 results show it is effective too. 

Here we have to say that this work is only about the estimation of un- 
known parameters and orders with the objective function (9) for for non- 
commensurate and hyper fractional order chaos systems in non-Lyapunov 
way. It can be concluded that DE in Algorithms 1 can be change to other 
artificial intelligence methods easily. For the cases that some fractional or- 
der differential equations are unknown but with definite orders q have been 
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discussed in Reference [36]. 

In the future, we will do further researches for the cases that neither the 
fractional orders nor some fractional order equations are known. That is, the 
objective function is chosen as (8) in the novel mathematic model in Section 
3. In another words, the objective will be changed into as equation (27). 

(?)/)* = argminF (27) 

(?,/) 

In conclusion, it has to be stated that proposed Algorithms 1 for fractional 
order chaos systems' identification in a non-Lyapunov way is a promising 
direction. 
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